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The blast wave resulting from an explosion was simulated to provide guidance for models 
estimating risks for human spacecraft flight. Simulations included effects of headwind on 
blast propogation, Blasts were modelled as an initial value problem with a uniform high 
energy sphere expanding into an ambient field. Both still air and cases with headwind were 
calculated. 


Nomenclature 


AEo Initial Explosion Energy Increase 
e Energy per Unit Mass 

y ^ Ratio of Specific Heats 

p Mass per Unit Volume 

Cp Specific Heats at Constant Pressure 

c v Specific Heats at Constant Volume 

m 0 Initial Explosion Mass 

p Pressure 

ro Initial Explosion Spherical Radius 

T Temperature 

Vo ‘^y SL Initial Explosion Volume 


I. Introduction 


Requirements for a abort system for a launch vehicle are driven in part by the necessity to escape a 
blast wave resulting from a system malfunction. Blast wave propogation is affected by the vehicle speed, 
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altitude, the extent of material involved in tne mast ana its spatial distribution. rreuiction oi blast wave 
overpressures and arrival times are important in deriving requirements for abort system thrust, trajectory 
and abort lead times. An investigation of blast wave propogation is ongoing as part of a probablistic risk 
assessment (PRA) project. 
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II. Method 


A. Euler Equations 

The Euler equations , written in conservation law form are 


Where 


9Q 9T\ 

— -t = 0 

9t 9xi 

(1) 

Q— [p, pui, pu 2 , pu 3 , pe T> ] 

(2) 

Fi = [pui, puf +p, puiu 2 , puiu 3 , pu - ) (e t + p)] 

(3) 

F 2 =[pu 2 , pu 2 Ui, pu|+p, pu 2 u 3 , pu 2 (£ T + p)] 

(4) 

F 3 =[pu 3 , pu 3 u,, pu 3 u 2 , puf +p, pu 3 (e T + p)] 

(5) 


where: £t = e +(uiUi)/2 



c v dT = c v T 


P = (Y-I)pe 


B. Numerical Method 

The code used in this study was OVERFLOW2. Matrix dissipation was used with smoothing parameters 
as recomended by earlier studies of high speed flows with this code. 1 The recommended eigenvalue limits 
(V £n = 0.3, Vet = 0-3) are not adjusted. 

Time accurate simulations were obtained with the dual time method. In general, simulations were run 
at the highest physical time step possible to allow the initial condition to be simulated stably. This is a 
conservative choice, as the simulation can be run at much higher time steps once the high variation initial 
transients integrated past. 

The relaxation method is the implicit Pulliam Chausee diagonal method, with timesteps chosen for a 
constant CFL (itime=3). CFL values for these high speed flows were chosen at 0.4, and the variable time 
step was adjusted as described in the overflow documentation. 2 


C. Scales, Nondimensional and Dimensional Coordinates 

A length scale of the problem is given by the cube root of the ratio of the explosion energy to the ambient 
pressure. Thus, A = (Eo/p a ) 1,/3 is a length which is a natural length scale of the problem. A natural velocity 
scale is given by the ratio of the initial energy to the initial mass. For example, a TNT explosion has a 
natural velocity scale of 2.044km/sec (4.18 x 1 0 6 ) 1 - /2 ). Pressures are non-dimensionalized by the ambient 
pressure. 

To translate a solution from the non-dimensional overflow coordinates, the length scale of the blast in 
physical coordinates is determined by the initial radius in computational coordinates and the corresponding 
physical extent of the initial sphere. Given the ratio of the explosion to ambient density, and the initial 
mass, the physical size of the sphere in meters is given by 


( Pa 3 nip \ 1/3 
\ Po 47rrho Q ) 


( 6 ) 


Thus, if the initial radius of the sphere is taken as unity in the non-dimensional computational coordinates, 
the physical length L corresponding with this length is given by this formula. The time scales (which in 
overflow are given by L/a.^ ) are then determined by this choice of L and the ambient sound speed. The 
reference Mach number is taken as unity, so that u® = a,*, . Pressures, temperatures and velocities follow 
the standard overflow methodology, where the ’infinity’ conditions are identified with ambient conditions. 


2 of 6 


American Institute of Aeronautics and Astronautics 



D. Iniitial Condition 


The initial condition was simulated as a sphere of high energy fluid at rest with respect to the inertial frame 
chosen for the simulation. Three parameters fully specify the conditions within the initial sphere containing 
the explosion, the additional energy AEo, the contained mass, mo, and the initial radius, To- Assuming a 
perfect gas, the total energy in a sphere is given by 


E 0 = 


47TT o t 
— — p 0 c v T 0 


— moc v To 


The energy(AEo) and mass(mo) are taken as known inputs to the simulation, and it is assumed that before 
the explosion the mass involved was already at ambient temperature, so that Eo = Po c vT Q 4- AEo, and thus 
the energy per unit mass (^- = e) is a constant for a given explosive yield and mass, If the perfect gas 
assumption is made, this impfies that the initial temperature of the ball is similarly a function of these two 
parameters. Given a more realistic assumption of chemically reacting gases, assuming chemical equilibrium 
and given element mass fractions, the temperature is similarly a function of these two parameters. 

The initial density of the ball is given by po — 47tmoro/3. given a fixed mass, then the density is a 
function of To , the initial extent of the explosion. The range of values plausible for density range from water 
on the high end(as kerosene is close to a specific gravity of one), down to the ambient atmospheric density. 
The high range corresponds to a condition where the fuel and oxidizer are perfectly mixed, and then ignited 
simultaneously. This is admittedly an unrealistic possibility, but is a plausible upper limit on the density 
likely to be encountered in any explosion. The lower limit would correspond to the fuel mixing and expanding 
before ignition to the ambient pressure (and having a molecular weight comparable to that of air), and then 
igniting simultaneously. Again, this is not a plausible scenario, but again provides a plausible lower bound. 
The sensitivity of the solutions to this range of initial densities is explored, and is found to be a relatively 
small source of variation in blast overpressure predictions. 

These blast calculations are intended to provide apply to conditions where a given mass of the fuel 
available is ignited, but not all the mass. The spatial extent of the fuel which is ignited then determines the 
amount of mass which should be used as mo which could range from a minimum of the amount necessary to 
provide the AEo (which assumes that all the fuel in a given volume is burnt along with the required amount 
of oxygen), to an essentially unbounded amount (allowing for the mixing of the fuel with the ambient air 
before ignition). The standard assumption used is that the total mass involved is equal to the mass of 
propellant available. 

Varying mo holding AEo fixed varies the temperature(which is given by ATo = AEo/(c v mo)). Holding the 
density fixed (so that r 0 varies) will then correspond to assuming that the density of the initial fuel/oxygen 
combination is roughly at fixed density before ignition, and varying mo is simply including more or less 
unburnt material in the initial explosion sphere. This type of parameter variation will be termed a ’constant 
density’ variation. 

Holding mo and AEo fixed, but varying ro corresponds to choosing varying initial densities, and results 
in a sphere with at fixed temperature (given by DeltaT 0 = AEo/(c v mo)) but with variation in the initial 
density p 0 — 3mo/(47tro). This parameter variation will be termed a ’constant temperature’ variation. 

Holding r 0 and AEo fixed, but varying mo results in a sphere at a fixed pressure (given by p = 3(y — 
1 )Eo/(47rr^)) but with varying temperature and pressure. This parameter variation will be termed a ’constant 
pressure’ variation. 


III. Results 


A. Comparison with Experiment 

One experiment which can be used for validation purposes for this case is the Boyer experiment. 1 This is an 
experiment with a glass sphere filled with either air or helium is shattered, This is a particularly valuable set 
of data because it has a well defined set of intial conditions. The pressure ratio and temperature ratio(which is 
unity) are known, and the initial condition is still air within the high energy sphere is zero velocity. A wealth of 

data is available for this case, including what are effectively experimentally determined x 1 diagrams which 

allow the blast wave speed predictions to be verified along with the overpressures as a function of distance. 
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The overpresssure for the case with air (7 = 1 .4) at a pressure 
ratio of 27.2(or 22 corrected) is compared with the results of the 
simulation in figure 1. The solid line is the line corresponding 
to the maximum pressure over the span of the simulation at the 
given spatial location. This was obtained by taking obtaining 
what would be equivalent to oscilliscope traces, the solution as a 
function of time at a given spatial location, then searching this 
trace for the maximum pressure. The experimental pressures are 
given by symbols, and snapsots of the pressure field predictions 
at given times are superimposed, showing the predictions for the 
blast wave propogation. The experimental results were corrected 
to account for the energy loss created by accelerating the shat- 
tered remains of the glass sphere. The experimental pressure 
measurements lie within the range of pressures bounded by the 
uncorrected and corrected pressure ratio results. 

The blast wave arrival time are another important prediction 
which the simulation will provide. The experimental time have 
a delay of 50 to 60 p seconds between the mallet striking the 
glass sphere and the emergence of the shock wave from in the schlieren record. It is assumed that this is 
the time required for the glass sphere to absorb the energy of the mallet and shatter. The arrival time for 
the blast wave at 5.7 inches is between 300 psec and 320 psec, which are predicted to be at 240 to 270 psec 
computationally. Given the delay time of 50-60 p seconds for the glass enclosure to shatter, this is good 
agreement with the experimental measurements. The arrival time at 8 inches is 440 psec experimentally, 
and this simulation prediction is between and 380 and 420 psec, again in good agreement when the glass 
shattering delay is taken into account. A possibly unwarranted extrapolation from the experimental x — t 
diagram has an arrival time at 10 inches of about 590 p seconds, which again is in good agreement, with the 
predicted 510 to 540 p seconds arrival time. 





Figure 1. Overpressure Predictions for Boyer Air 
Blast 


(a) Primary Blast Wave Arrival 


(b) Primary /Secondary Blast Wave Arrivals 


Figure 2. Comparison of Arrival Time Predictions with Boyer Experiment 

Another ’arrival time’ prediction is for the time delay between the arrival of the primary and secondary 
shock at the r= 5 . 7 in postion. The time delay between the arrival of these two shocks may be obtained the the 
experimental oscilliscope trace(Figure 2 (a)), and is 580 psec experimentally. The computational oscilliscope 
trace is shown in figure 2(b), and the delay between 580 psec and 610 psec. 
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B. Explosions in Still Air 

The initial solutions obtained were obtained for still air conditions. For an explosion with an initial energy 
of 1KT of TNT (equivalent to 3.8 terajoules), with a mass of 5KT(note that the mass is assumed to be 
five times the mass of the TNT equivalent), this yields an internal energy of 836 KJoules/kg. Assuming a 
molecular wieght of 29 kg/kMole, a ratio of specific heats of 1.4, this yields a temperature increase of 


AT 


AE 0 (Y-l )M 


m 0 
836 x 10* 
11 64.3° C 


n 

29(0.4) 


8314.3 


°C 


again assuming perfect gas conditions. Note that lowering y or increasing the mass decreases this temper- 
ature increase. For instance, if the same calculation were done for TNT, where the mass was 1KT, but y 
was decreased to 1.3(to account for non-perfect gas behaviour), the corresponding temperature increase is 
4,545°C, substantially hotter. 

The overpressure for the 1KT TNT energy with 5 KT mass case in universal coordinates is given in figure 
3(c) for a variety of density assumptions. The resulting variation in overpressure in confined to the nearfield, 
and is much smaller than the corresponding overpressure variation shown by a variation in either explosive 
yield or mass in figure 3(b). 





(a) Constant Initial (b) Constant Initial (c) Constant Initial 

Pressure(Eo , tq Constant) Density (mo , tq Constant) Temperature(Eo , mo Constant) 


Figure 3. Overpressure Predictions for nominal 1KT TNT Explosion, with 5KT Mass 


C. Explosions in Headwinds 

The impetus for this study was to provide an check on estimates of the effects of headwind on blast pressures. 
Representative plots of overpressure predicitions for blasts occuring in the presence of a headwind are shown 
in figure 4. Figure 4(a) shows the effect of lowering the mass of the initial blast, keeping the explosive 
energy constant. This yields higher overpressures initially, but the shock wave actually does not extend as 
far upstream as the case with the higher mass. Figure 4(b) shows the more subtle effect of varying the 
spatial extent of the initial blast, holding both blast energy and initial mass constant. There are some initial 
differences in overpressure, but they are much less pronounced and both cases ’turn around’ at the same 
location. 

One of the features of the case with headwind not present in the still air case is the finite extent of 
influence in the upwind direction of the blast. This critical difference is one of the reasons for doing this 
study, which is to better elucidate the escape system requirements for potential launch vehicles. 


5 of 6 


American Institute of Aeronautics and Astronautics 






r 


556MJ Explosion at 11.58km Altitude, M = 1.5 Headwind, 6.4m Radius 



11.58km Altitude, M=1.5 Headwind, 556MJoules, 1.1x10® kg 



(a) Varying Initial Mass 


(b) Varying Initial Density 


Figure 4. Overpressure Predictions for Nominal 141 Ton TNT Explosion with Headwind 


IV. Plans 

Reacting flow simulations are planned to elucidate some rudimentary effects of finite rate chemistry and 
caloric imperfection on the blast wave propogation. Conditions at various altitudes, various headwinds, 
and various explosive yield/mass combinations are being investigated, and will be reported. Overpressures 
produced on vehicles will be investigated to further elucidate the blast wave effects on crew safety. 
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